On the Asymptotic Formula for the Number of Plane 
Partitions of Positive Integers 

Ljuben Mutafchiev 

American University in Bulgaria and 
Institute of Mathematics and Informatics of the Bulgarian Academy of Sciences 

1 j uben@aubg . bg 
and 

Emil Kamenov 

Sofia University 
kamenovQf mi . uni-sof ia . bg 

A plane partition uj of the positive integer n is an array of non-negative integers 

(1) 



1^1,1 ^1,2 t^l,3 
^2A 1^2,2 1^2,3 



for which 



and the rows and columns are arranged in decreasing order: 

for all i,j > 1. The non-zero entries ujij > are called parts of uj. Sometimes, for the sake of 
brevity, the zeroes in the array (1) are deleted. For instance, the abbreviation 

3 2 1 
1 1 

is assumed to present a plane partition of n = 8 having 2 rows and 5 parts. 

It seems that MacMahon was the first who introduced the idea of a plane partition; see 
[5]. He deals with the general problem of such partitions enumerating them by the size of 
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each part, number of rows and number of columns. These problems have been subsequently 
reconsidered by other authors who have developed methods, entirely different from those of 
MacMahon. For important references and more details in this direction we refer the reader to 
the monograph of Andrews [1; Chap. 11] and to the survey paper of Stanley [8; Chap. V]. 

Let q{n) denote the total number of plane partitions of the integer n > 1. It turns out 
that 

oo oo 

Qix) = 1 + E ^H^" = 11(1 - x^)'^ (2) 

n=l j=l 

(see [1; Corollary 11.3] or [8; Corollary 18.2]). The asymptotic of q{n) has been obtained by 
Wright [10]. It is given by the following formula: 



:[7o + E7.n-^'/' + 0(n-2(^+^)/=^)], (3) 

where 



XI, „ . 

h=l 



is the Riemann's zeta function. 



y^ogy 

dy (4) 



Q2ny _ I 

and 7/i are certain constants with explicitly given values. © implies the asymptotic equiva- 
lence 

^(-) - lH!^""'^'' {3[C(3)] V3(,/2)2/3 + 2c}, n oo. (5) 

Many authors refer to (jH)) and (0) in their work on this subject. Wright [10; p. 179] states 
in his theorem that 70 = 1, however, at the the end of its proof his argument shows that the 
true value of 70 is 70 = 3^^/^. In fact, the substitution v = 1 + y [10; p. 188] that he makes 
in the integrals denoted by generates the factor 

(3 + 2?/)'(™+^/2)^^^ 0,1,... (6) 

in their integrands. Then, Wright shows how determine the coefficients 7^ in Q. He 
proves that = (~l)"^^sm/2'/r, where bsm is the coefficient of y"^"^ in the expansion of 

(l+^)2«+2-+13/12(3^2t/)-™-l/2 

(here s is a non-negative integer; see also the definition of the constants on p. 179 [10; 
formula (2.15)]). It turns out that ii h = 0, one has to set in © m = 0, and thus by [10; 
formula (2.22)], 70 equals the constant term in the power series expansion of 

(l + yfs+i3/i2^^^2y)-'/\ 

Hence 70 = 3'^/'^. 
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The asymptotic equivalence (0) can be also obtained after a direct application of a theorem 
due to Meinardus [6] (see also [1; Chap. 6]) who has obtained the asymptotic behavior of the 
coefficients in the power series expansion of infinite products of the following form: 

oo 

where {cij}j>i is given sequence of non-negative numbers. He introduced a scheme of as- 
sumptions on {aj}j>i, which are satisfied by the generating function ^ of the numbers g(n). 
Further generalizations of Meinardus' result are given in [3]. 

The object of this work is to show that Meinardus' theorem implies asymptotic formula 
© with 70 = 3^1/2. This fact was also briefiy mentioned in [4], where the limiting distribution 
of the trace of a random plane partition was studied. A sketch of the proof of (jSI) based on 
the fact that Q{x) is an admissible function in the sense of Hayman [2] was presented there. 

We start with the statement of Meinardus theorem [6] (see also [1; Chap. 6]) . Let 

oo _ oo 

n (l-x^)""' = 1 + 5]r(n)x", 

j=l n=l 

where x = , 3?e r > 0. Consider also the associated Dirichlet series 

°° a ■ 

^(^) = E^> s = cT + tt. (7) 

j=i J 

Meinardus assumes that the following four conditions hold. 

(i) D{s) converges in the half-plane a > a > 0. 

(ii) D{s) can be analytically continued in the region 3fJe s > —Co (0 < Cq < 1), where 
D{s) is analytic except for a simple pole at s = a with residue A. 

(ill) There exists a constant Ci > such that D{s) = 0(|t|*"^) uniformly in a > — Cq as 
\t\ oo. 
(iv) Let 

oo 

fi'(^) = Z!%^^ x = e"^ (8) 

where v = y + 27riw and y and w are real numbers. For | argf | > ir/A and \w\ < 1/2 and for 
sufficiently small y, g{v) satisfies 

^eg{v)~g{y)<-C2y-', (9) 

where e > is an arbitrary number, and C2 > is suitably chosen and may depend on e. 
Theorem [6] If conditions (i) - (iv) hold, then 

r(n) = Cn^exp|n^ + ^) [^^(a + l)C(a + (1 + ©(n"^^)), 

where T{a) is Euler's gamma function, 

C = e^'(°)[27r(l + a)\-^'^[AT{a + l)C(a + 1)]^^, (10) 



3 



and 6 is an arbitrary positive number. 

In the case of plane partitions we have {aj}j>i = {j}j>i and thus by ((7j), D{s) = ({s — 1). 
Thus, (i) - (iii) follow from well known properties of the zeta function (see e.g. [9; Chap. 13]). 
Moreover, in (ii) Co is an arbitrary constant within range (0, 1), since the zeta function has 
analytical continuation in whole complex plane. It is easily seen that a = 2, A = 1. To verify 
condition (iv) notice first that for {aj}j>i = {j}j>i the function g{v) defined by (jHI) becomes 

oo 

i=i 

Hence 

^egiv)-giy) = ^ je-^^cos (27rj^) - 1] = T.i-'^) -77^T.f''''^'''■ 



j=l k=l V-'V j=l 

We simplify the inner sum here replacing its terms by the following derivatives: 

2k+l -jy _ (_i \2k+l_^ (p-m 



Further, to calculate the {2k + l)th derivative of l/(e^ — 1) we apply Faa di Bruno's formula 
[7; Section 2.8]. Using also an identity for the Stirling numbers of the second kind o"^™"^ [7; 
Section 4.5], we find that: 

d^'^+i / 1 (-l)"^m! ^ {2k + l)\e"'y 



2k+l ( '\\m^\pmy 2k+l 

m=0 \^ ^) m=0 



where J2b means that the summation is over all non-negative integers bj that satisfy n = 
J2j jbj, m = J2j bj- It is clear that all terms in the last sum are close to zero if y is sufficiently 

small except the last one. It equals {-iy^+\2k + = (-l)^*^+^(2fc + 1)!. Therefore, 

after simple algebraic manipulations, we can rewrite (fT^ in the following form: 

oo 

^egiv) - g{y) = Y.{-l)\2'nwf\ey - l)-^^^^^^[2k + 1 + Mv)] 

k=l 

■[3 + ^i(2/)]+ -^-^[5 + ^2(2/)]-7^— ^[7 + V>3(?/)] +..., (14) 
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where 

My) = E i-irm\e-y{ey - ir^'--a^£, = 0{y) 

l^/^J- m=0 

as y ^ O"*". For sufficiently small y we also have e^ — 1 = y + 0{y'^). Moreover, the requirement 
I argr |> 7r/4 implies that > y^^. Using these arguments, we conclude that the first 

term in (14) is — (3/?/^)[l + 0{y)] as y — > O"^. For the other terms in the curly brackets of (jl4j) 
we get the estimate 

+ 1 - + 3 + 0(,)|} < ^|-2 + 0(,)l < -2,- + Ofe-). (15) 

where / = 1, 2, .... Suppose now that e G (0, 2]. Estimate (15) implies that one can find certain 
constant C2 such that 

y^[^egiv)-giy)]<-C2<0. 

If e > 2, then for any rj > and sufficiently small y we have < y^"^ < rj. Therefore for such 
y's one can define a sequence of enough large positive integers {My} satisfying the inequality 
rj/2My < ?/^~^, or equivalently, —77 > —2Myy''~^. Now, estimate (fT3j) shows that —rj exceeds 
the sum of My terms of the expansion of y'^['Reg{v) — g{y)] and this sum in turn exceeds the 
sum of all terms of the expansion. Hence, for e > 2, we also obtain inequality Q of condition 
(iv) with C2 = Tj. 

To calculate the constants -D(O) and -D'(O) we use the following functional equations (see 
[9; Chap. 13]): 

C(l -z)=2 cos ^TTZ {2'r)-^T{z) C{z), (16) 
r(^) C(^) = / -z-^dw = {2ny / -^—-dw. (17) 



e"^ - i Jo e^™ 



Substituting 2; = 2 in (|TI)|) we get: 



D(o) = c(-i) = -^. (1^ 



Next, differentiating (fTB|) and (fT7j) and combining the expressions of the corresponding deriva- 
tives, we obtain 

aw. 



C'(l - z) = sin -nz 2"%^^^ r(^) C(^) - 2 cos -tiz f 
2 2 Jo 

Substituting again z = 2, we find that 



D'{0)=a-l) = 2l ^^^dw = 2c, (19) 



where c is the constant defined by (4). 

From (Uni) - (HI), (HHI) and (Uni) we obtain 
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7/36 „2c 



K 



36' 3 

where ei > 0. Notice that in (fT^ we assume that 6 is sufficiently small and Cq is close to 1. 
Then, applying Meinardus theorem, we obtain 



q[n) 



icm 



7/36 



n 



-25/36 



exp {3[C(3)]i/=^(r2/2)2/3 + 2c} 



211/3631/2^1/2 

which shows in turn that Wright's formula [10] with 70 = 3"^^^ is valid. 

We also present some numerical computations of g(n) based on formula (jSJ and the fol- 
lowing recurrence: 

n 

n q{n) = ^ q{n — k) /?2(^)- 



k=l 



Here P2{k) is the sum of the squares of the positive divisors of k (see [1; Section 14.6]). We 
used Maple to get the results summarized in the table below. The exact values of q{n) are 
listed in its second column. The third and fourth columns contain the corresponding values 
of q{n) computed when 70 = 1 and 70 = 3^^^^, respectively. 



n 


q{n) 


7o = 1 


70 = 


10 


500 


910.69 


525.79 


100 


59 206 X 10^2 


103 709 X 10^2 


59 876 X 10^2 


1000 


35 426 X 10*° 


61 507 X 10*^° 


35 511 X 10*° 


10 000 


45 075 X 10^^^ 


78 113 X 10^^^ 


45 098 X 10^°^ 
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